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Abstract 

A new result for the hadronic vacuum polarization correction to the muonium 
hyperfine splitting (HFS) is presented: 

Ai/(had - vp) = (0.233 ± 0.003) kHz. 

Compared with previous calculations, the accuracy is improved by using the latest 
data on e~^e~ -^ hadrons. The status of the QED prediction for HFS is discussed. 



1 Introduction 

Hyperfine structure of two-body atomic systems is of interest for a variety of reasons. 
In particular, it offers opportunities for precise tests of bound-state QED and accurate 
determination of fundamental constants like the fine structure constant a, or electron- 
to-muon mass ratio. The hyperfine interval in the hydrogen ground state used to be 
the most precisely measured quantity. The fractional uncertainty was about 10~^^ (see 
e.g. DH), although the theory cannot make a prediction with accuracy better than 10~^. 
The theoretical uncertainty comes from the proton magnetic form factor unknown at low 
momentum transfers [0 and from the proton polarizability contribution. 

Studying purely leptonic systems such as muonium (the bound state of a positive 
muon and an electron), one can avoid problems of the proton structure. Despite the short 
muon lifetime (about 2.2 yusec), the hyperfine splitting (HFS) of the muonium ground 
state has been measured very precisely |0, ^ , 

z/HFs(exp) = 4463 302.776(51) kHz, (1) 

with a fractional uncertainty of 0.011 ppm. 

The theoretical prediction for the HFS is determined by the Fermi energy, resulting 
from a non-relativistic interaction of electron and muon magnetic moments. It can be 
expressed as a combination of fundamental constants 



Up = — (^«) cZ Hoc — 



1 + a^) . (2) 



m 
Another convenient formula which will be used to determine the Fermi energy is 



vp = -7r{Za) cZ R^—- 
3 ^1% 
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(3) 



The accuracy of these expressions is limited by the electron-to-muon mass ratio m/M or 
the ratio of the muon magnetic moment to the electron Bohr magneton /i^//i^. Since 
the muon anomalous magnetic moment a^ is known with high accuracy [^ |^, 0]) these 
two formulae are essentially equivalent. Given that the three most precise experiments 
0, H, §] determined the magnetic moment ratio, we will use the equation (|^). Using 
the values 1/a = 137.035 999 58(52) 0, cRoc = 3 289 841960 368(25) kHz 0, /i^//ip = 
3.183 345 24(37) [|, |, and ^p/fx% = 1.521032 203(15) ■ lO"'^ g, we find 

Up = 4459 031.920(511)(34) kHz, (4) 

where the first (larger) uncertainty is due to the magnetic moment ratio and the second 
to the fine structure constant. The value M/m = 206.768 276(24) used in this work was 
obtained by us from the values of the magnetic moments above as well as the experimental 
value a^ = 1.165 920 3(15) ■ 10~^ [^. Here Z is the nuclear charge and in the case of 
muonium Z = 1, however, it is convenient to keep that for the classification of different 
contributions (see Section H). In practice, precise measurements of the muonium HFS 
P, M serve as the most precise determination of the muon-to-electron mass ratio P| . In 
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Figure 1: Hadronic contribution to muonium HFS via the muon anomalous magnetic 
moment. 





Figure 2: Bound-state hadronic effects to muonium HFS considered in this paper. 



the future, muon mass might be determined more precisely from other experiments |T0| 
thus enabling a new test of bound-state QED. 



Further experimental progress in the muonium HFS is also expected |]T0[ if high- 
intensity muon sources become available. In view of those advances, it is conceivable that 
the accuracy of bound-state QED tests may become limited by interactions beyond QED, 
in particular by the strong interaction effects. In this paper we study hadronic vacuum 
polarization, which is the leading contribution of the strong interactions. These effects 
influence the HFS in two ways: through effects included in the Fermi energy (contributions 
to the muon a^ shown in Fig. p and bound-state effects depicted in Fig. ^. Hadronic 
effects in a^ are relatively large (~ a^) and have been studied by a number of authors 
over the last 3 decades (see references in P, |ll|). In the present case, we are using an 
experimental determination of the muon magnetic moment (/i^) and do not address this 
theoretical issue. We focus on the bound-state effects of the hadronic vacuum polarization 
(~ a^{m/M)). 

The theoretical expression for the hyperfine splitting can be written in the form 

z/HFs(theor) = vp + Az/(QED) + Az/(weak) + Az/(had). (5) 

The biggest correction to the Fermi splitting {vp) comes from the QED effects (Az/(QED)) 
and it is under consideration in Section ^ The two terms beyond QED are due to weak 
(Az/(weak)) and strong (Az/(had)) interactions. 

Influence of the weak interaction effects on the muonium HFS was studied in the 



leading order |T^, |13|. The leading correction is induced by the neutral currents, given by 
a Z-boson exchange diagram. Recently it was found that the sign of the corrections had 
not been well understood and the absolute value was verified. In a recent paper |T^ that 



problem was solved. The contribution is 

« / N Gp-mM 3 
Az/ weak = ^ —- ■ up ~ -0.065 kHz , 6 

where Gp is the Fermi constant of the weak interaction. Next-to-leading-order contribu- 
tions, studied in Ref. [|T3|, are O (1%) of the leading contribution (^ and thus negligible. 



The strong interaction effects (Fig. |]) were studied in [|l^, |T^, Q. We examine the 



hadronic contribution in Section ^following the approach of Ref. [0], in which the QED 
and hadronic parts are separated and the final result in the leading order is presented in 
the form 

AKhad - vp) = - — - vp\^^^^ d. a{s) H{s), (7) 

where s is the center-of-mass energy squared, a{s) is the total cross section of electron- 
positron annihilation into hadrons and H[s) is the QED kernel calculated in ^^. Our 
calculation gives 

Ai/(had - vp) = (0.233 ± 0.003) kHz. (8) 

The uncertainty of this result, based on recent electron-positron annihilation data, is less 
than half of that of the most accurate previous calculation [|l^. The higher-order hadronic 



corrections are considered in part in |]T9 
In Section 13 we discuss our results. 



2 QED contributions 

The problem of the theoretical accuracy of QED effects in the muonium HFS was con- 
sidered in [|15| , where the importance of higher-order logarithmic corrections was pointed 
out. The complete QED theory includes a variety of contributions: electron and muon 
magnetic moments, external field terms, recoil and radiative-recoil corrections. The theo- 
retical expression can be presented as an expansion in small parameters: a, Za and m/M. 
In the case of muonium Z = 1, however, it is customary to keep Z in order to distinguish 
between QED effects (which are counted by powers of a) and the binding Coulomb effect 
which leads to an appearance of Za. In some papers the muon charge is denoted by Ze 
and so the QED effects due to the radiative corrections on the muon line involve Z'^a. 
Since that is still a QED effect, in our paper this correction is just a. 

Here we discuss briefly the progress made since the pubhcation of Ref. [jT3|. We divide 
the QED effects in three parts, 

Az/(QED) = Az/(ae) + Az/(QED3) + Az/(QED4), (9) 

because of their different status. The first term, 

Az/(ae) = aeUF = 5 170.926(1) kHz, (10) 



arises due to the electron anomalous magnetic moment Oe = 1 159 652 187(4) ■ 10 ^ for 



which we use the experimental value pO[. The next term, z/(QED3) includes the bound- 
state QED effects up to the 3rd order in the relevant small parameters, to be discussed 
in Section p.l| . These effects are known with good accuracy. The last term, z/(QED4), is 
related to the 4th order contributions and is subject to ongoing investigations. Its present 
status will be reviewed in Section 



2.1 QED up to third order 

The contributions up to the third order in those parameters give 

Az/(QED3) = 

= upl-iZaf + a{Za)(\n2- 



a{ZaY 
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Origins of individual terms in Eq. (|TT]) are discussed in detail in the review 
numerical value resulting from ([TT|) is 

Az/(QED3) = -899.557 kHz, 

with an error beyond the displayed digits. 




(12) 



2.2 Fourth order 

The essential part of the fourth order corrections has been evaluated (see e.g. ||2^). The 
known terms are summarized in Table |I[ Some of them were calculated in the leading 
logarithmic approximation (using In 1/a ~ In M/m ~ 5) and in those cases we follow [|1^ 
and estimate the uncertainty due to the non-leading terms as half of the leading ones. 
Table ffl contains a number of recent results which appeared after the review [|15l. Some 



of the contributions will be reviewed in more detail below. We also present in that Table 
a very small contribution of the tau lepton loops replacing the hadronic loop in Fig. |^. 
This effect is O {a{Za)mM / rn^) and amounts to about 2 ■ 10~^ kHz. 



2.3 One— loop corrections 

The electron one-loop self-energy has been studied in some detail. After direct calcula- 



tions of the a{Za) term |3^, ^, two independent exact (without an expansion in (Za)) 



Correction Contribution to the HFS 



Reference 



{Za) 



Total 



0.03 kHz 



{Zaff: 


-0.29(13) kHz 


(W(f)' 


-0.02(1) kHz 


(^«)(f)' 


-0.02 kHz 


a{Zaf 


-0.52(3) kHz 


^{Zaf^ 


0.39(17) kHz 


<Za){^)' 


-0.04 kHz 


a\Zaf 


-0.04(2) kHz 


o^\Za)f^ 


-0.04(3) kHz 


a\Za) 


±0.01 kHz 


T lepton 


0.002 kHz 



-0.55(22) kHz 



mm 



m 



see Table I 



24, 26, 27 



m.m 



1 24 



see caption 



Table 1: The fourth order corrections. The first column gives the order of corrections 
relative to up. The uncertainty due to the unknown term a^{Za) is estimated as (a/vr)^ 
times the known a{Za)i'F term in eq. (pHI). 



numerical evaluations were performed [^^ ^ . Unfortunately, the accuracy for Z = 1 was 
not high in both cases and the final results for muonium were obtained by extrapolating 
the data from higher values oi Z to Z = 1, using previously found coefficients of a{Zay 
and a{Za)^ In Za []T3[. The result of the later calculation |5^ is presented in Table 



1331 



^ as a{Za)^^. It contains the non- logarithmic a{Za)^ term and higher-order corrections. 
The result of [^ is slightly higher than that of [^. The calculation gives for this term 
■14:.3{1.1) alZaYup/iT [Q and —12{2)a{Za)^i'F/'7r |^^. Recently the non-logaricthmic 



term was calculated directly by expansion in Za to be —15.9{1.6)a{Za)^i'F/7T [40 

The uncertainty of the extrapolation [Q is determined as an estimate of the unknown 
higher-order terms rather than the accuracy of numerical data. However, the final un- 
certainty is not clear, because of the value of the non-logarithmic a(Za)^ term p5| , ^ 

(cf. m my- 



:i7.122 339...)-^^^^i^^ [35] 



TT 



TT 



, aiZa)"^ , — K 

:i7.122 7±0.0011)- ^ -Up !§ 



(13) 



These results appear to be in fair agreement with each other. However, it has been 



pointed out [^ that the term-by-term comparison reveals inconsistencies. The low- 



energy contribution in Ref. ||36| is by 0.0257(2) higher than that in Ref. |^, while the 
medium- and high-energy parts are by 0.0246(14) lower. It should also be mentioned that 
different fits for extrapolation to the zero photon mass A in Ref. [^ are not consistent. 
Such extrapolation was necessary to calculate the high momentum part. These problems 
remain to be clarified. Since the same method was used in [^], it is the result of [^ that 
we include to Table 0. 



Order 



Contribution to the HFS Reference 



Comments 



a{Zaf In^ 
a{Zay+ 

a{Zaf In^ 
a{Zaf 
a{ZQif 



-0.53 kHz 

-0.06 kHz 

0.03 kHz 

0.03 kHz 

±0.03 kHz 



in 



ii 

m 
m 



see Sect. |2 



self-energy 

self-energy 

vacuum polarization 

vacuum polarization 

Wichmann-Kroll 



one- loop 



-0.52(3) kHz 



sum of all terms 



Table 2: Higher-order one-loop contributions. 



Recently, the one-loop vacuum polarization was calculated exactly (to all orders in 
Za) W\ B3. For our purpose it is enough to know the aiZaf terms only 



13 



539 \ 



l\v = aiZaf — In \ ] up 

^ ' V24 Za 288/ 



(14) 



the logarithmic part of which was obtained in Ref. [T^, while the constant was found in 
Ref. [||. 

The only unknown term in the order a{ZaY is now the so called Wichmann-Kroll 
contribution. We estimate it by the value of the non-logarithmic part of the VP term in 
Eq. (inp. 



2.4 Recoil effects 

A number of results on various recoil effects were obtained after the publication of |l^ . 



Pure recoil effects of the order {ZaY were studied without an m/M expansion in 
Ref. [Q for arbitrary mass and numerical results for several values of m/M were 
obtained. Recently the same correction was found for muonium [^ (see Table |ip. 

For radiative recoil corrections {a{Za)m/M), there are minor discrepancies between 
the published numerical and analytical results, summarized in Table ^. All analyt- 
ical results were obtained by two groups which are in agreement. It is likely that 



the uncertainty of the numerical calculations was underestimated [Q] and the dis- 
crepancies have no connection with higher-order corrections. This explanation is 
supported by the similar situation with the radiative recoil corrections in positron- 
ium, shown in the same table. 



Contribution 



Numerical 
result 



Analytical 
result 



Discrepancy 



e— line 

/i— line 

Positronium 



3.335(58) H, |T6| 
-1.0372(91) 
-1.787(4) 



^m 



^m 



3.499 m, 13, |T], |2 
-1.0442 [|8],|1|] 



1.805 H, Ig, |5T 



0.164(58) 
-0.0070(91) 
-0.018(4) 



Table 3: Radiative-recoil corrections in units of a^{m/M)vF- The numerical results 
include all orders in m/M. The analytical results refer to sums of the third and fourth 
order, a'^[m/M) and a'^{m/M)'^ for muonium. For positronium, M = m and up is defined 
in m. 



The uncertainty of the QED calculations of the muonium hyperfine structure arises 
mainly from two sources: a{ZaYm/M and {Za)^m/M. Only double logarith- 
mic corrections were calculated until recently, when some single-logarithmic terms 
(~ ln{Za)) were found. The single-logarithmic terms with a recoil logarithm 
(ln(M/m)) are still unknown. We discuss a possible consequence for the accuracy 
of the theoretical calculations in the next section. 



2.5 Non-leading logarithmic terms and uncertainty of the cal- 
culations 



In this paper we follow |T5| and estimate any unknown non-leading terms by a half-value 
of the corresponding leading logarithmic terms and sum them as independent. There 
is no proof that the non-leading terms can be estimated in this way. However, all ex- 
perience with different contributions to the Lamb shift and the hyperfine structure (see 
e.g. the review |^| with a collection of various terms) demonstrates that this assumption 
is reasonable. 

The fourth order leading corrections contain cubic and quadratic logarithmic terms 
and in the case of some contributions, known in the logarithmic approximation, more 
than the leading term has been evaluated (see Table |1| for references): 



a'^{Za) m 



-jj-O 



M 



Up X 



4 o M 4,2^ 

-- In^ h - In^ — 

3 m 3 m 



aiZa)^ m 



TT 



M 



Up X 



16 , 2 1 / 32 , 

— n^ + 1 + — n 2 

3 Za \ 3 



90/ Za} 



and 

^ ^ Up X -3 n— n n^ + 201n2 n \ . 

71 M \ m Za 3 Za \ 9 J Za) 

Even in this case we estimate the unknown terms as the half-value of the leading term. 
From comparison with e.g. known contributions of the third order a(Za)^z/i? (see Eq. 
(p!T|)) one sees that the constant term can be even bigger than the non- leading logarithm. 
Error estimates based on a non-leading logarithmic term may be misleading. The reason 
is that the leading logarithm in most cases originates from relatively simple diagrams 
under well understood conditions. The cancelations between different contributions occur 
very seldom. The leading logarithmic term has a natural value and is useful for estimates. 
In case of a cancelation which, although rare, is possible, an estimate should be based 
on the contributions before the cancelation. In case of the non-leading logarithmic terms, 
there are usually a few sources of different nature and some cancelations take place very 
often. 

The total magnitude of the muonium HFS interval calculated within QED 

i/QED = 4463 302.738(511) (34) (220) kHz. (15) 

is found by adding the values given by Eqs. (|10[), ([T^), and the sum of the 4th order 
contributions listed in Table |I]to the Fermi splitting (H). The third uncertainty is due to 
the 4th order QED effects in Table ^ 

3 Hadronic contributions 

3.1 Calculation of the hadronic contributions to muonium HFS 

The lowest order hadronic contribution to the muonium hyperfine splitting is given by 
the following expression: 

1 Tfl f°° 

Az/(had - vp) = -^jT ^F ds a{s) H{s), (16) 

ZVT'' M J Ami 

where the kernel His] calculated in [O is: 



rr/ X f s ^\ , ^ + r f s 3\, s 1 / 4M2 ^^_ 

His) = r + 2 rn - + - n — ^ + -, r = \l .(17) 

^' V4M2 J 1-r V4M2 2/ M^ 2' V s ^ ' 

The bulk of those effects is computed using the experimental data on the cross sections 
of e+e~ — i> hadrons in the energy range ^/s < y^ where sq is a scale above which pertur- 
bative formulae can be used. We choose ^/so = 12 GeV. After that one performs direct 
numerical integration of the experimental points similar to the approach of [^ where 
hadronic corrections to a^ were calculated. In contrast to the methods in which some ap- 
proximation of the data is used for the integration, in our approach model dependence is 
avoided as much as possible and, moreover, the calculation of the uncertainties is straight- 



forward. In addition to the data set used in [^, one can take into account significant 

9 



progress in the measurement of the hadronic cross sections in the energy range below 1.4 
GeV with two detectors at VEPP-2M in Novosibirsk [p3| , [5^ and the R determination in 
the energy range 2 to 5 GeV by the BES detector in Beijing [^ . 
The integration procedure gives for the contribution of this part 

Az/(bulk) = (0.2031 ± 0.0031) kHz. (18) 

The narrow resonances {u, 0, J/ip- and T-famihes) are evaluated separately. In a zero 
width approximation the contribution of a resonance B with mass rriB and electronic 
width Fee is 



Az/(resonance B) 



6 m T{B 



e ' e 



rriB 



-H{ml) 



Aa{m 



A"W - ^i:«}iv,(,„-i,4). 



(19) 



The leptonic widths relevant for our study should correspond to the lowest order 
(Born) graphs. However, experimentally measured leptonic widths listed and averaged 
in the Review of Particle Physics [Q, contain an additional contribution of the vacuum 
polarization by leptons and hadrons. Therefore, for the transition to the lowest order 
widths one should multiply the experimental values by [1 — Aa(mg)] (see the discussion 



of this issue in [^). Qf and ruf are the charge and mass of the fermion /, and N^f 
is the number of colours for the corresponding fermion {Ncf = 1 for leptons). The 
above approximate formula for Aa{s) is valid for s <C Mw and describes contributions of 
fermions much lighter than J s/A. In the present calculation we consider only effects of 
the lightest leptons, e and /i, and neglect the hadronic contributions to Aa(s). Hadronic 
loops are relevant only for heavy flavours (the J /ip- and T-family of resonances), which 
however give a smaller contribution and have a larger relative error. The hadronic effects, 
had they been included, would have shifted the J/?/^-family contribution by a few percent. 
The T-family contributes about 2 ■ 10~^ of the total resonance contribution. 
The sum of the individual resonance contributions gives 

Az/(res) = (0.0290 ± 0.0006) kHz. (20) 

In the region above sq one can use a perturbative formula for the hadronic cross section 

A-na^ 



(Tis] 



3s 

R{s) = R^^\s 



■R{s), 



as „ /«s 
1 + — + C2 — 

IT \ TT 



IT 



R^'\s) = ^c QiVi/3 + ^iV^/a) = y is>So) 



f2r 



where A^i/3(2/3) is the number of "active" quarks with a charge 1/3 (2/3) at given s, 
C2 = 1.411 and C3 = -12.8 E7[ 



At s ^ M^ one can use the asymptotic formula for H{s) [T^: 
H(s)- 



M2 /9 , s 15\ 



f22^ 



10 



so that after the numerical integration one finally obtains 

Az/(cont) ~ 0.0012 kHz. 



(23) 



The final result for the lowest order hadronic contribution is found as a sum of 



eqs. (^j 

Az/(had - vp) = Az/(bulk) + Az/(res) + Az/(cont) = (0.233 ± 0.003) kHz. (24) 

The individual contributions are presented in Table ^. It is clear that the dominant 
contribution (about 84%) comes from the very low energy range below 1.4 GeV studied 
in Novosibirsk. The regions below 1.4 GeV and from 1.4 to 3.1 GeV give about the same 
contribution to the total uncertainty. For most of the hadronic channels below 1.4 GeV 
analysis is still in progress in Novosibirsk. However, one can hardly expect significant 
improvement in this region since the dominant fraction of the error comes from the 27r 
channel which cross section is already known with a very high accuracy of about 0.6% 
around the p meson |^ . Very promising for decreasing the uncertainty of the 27r channel 
could be use of the r lepton data [^, but the real accuracy of this approach has been 
recently questioned |^ ^. The improvement of the uncertainty for ^/s from 1.4 to 3.1 
GeV will require new e"'"e~ colliders covering this energy range [^. Some progress in this 
region is also possible due to the B-factories and CESR which can contribute by using 
hadronic events from initial state radiation |]63| , |64 |. 



Final state Energy range, GeV Az/, Hz 



27r 


0.28-1.4 


158.8 ±1.9 


LO 




12.4 ±0.4 







13.2 ±0.4 


Hadrons 


0.6-1.4 


10.7 ±0.8 


Hadrons 


1.4-3.1 


23.8 ±2.2 


J/^ 




3.4 ±0.2 


Hadrons 


3.1-12.0 


9.8 ±0.5 


Hadrons 


> 12.0 


1.2 



Total 



233.3 ±3.1 



Table 4: Contributions to muonium HFS 



3.2 Comparison with other calculations 

It is instructive to compare the relative contributions of different energy ranges for i^hfs 
and a^. From Table | one can see that they are very close to each other indicating 



11 



the importance of the low energy regions. This is a natural consequence of the similar 
kernel structure, so that the main contribution in both cases comes from the low range 
of y/s. One can also note some enhancement of the high energy contribution to the 
HFS compared to that to the anomalous magnetic moment. That is due to the different 
asymptotic behaviour of the QED kernels for these two problems. For a^ the asymptotics 
of the kernel is proportional to M^/s, whereas for the muonium HFS it contains an 
additional logarithmic enhancement. 



Final state 


Energy range, GeV 


Az/, % 


O/,, % 


27r 


0.28-1.4 


68.0 


71.8 


UJ 




5.3 


5.7 







5.7 


5.8 


Hadrons 


0.6-1.4 


4.6 


4.2 


Hadrons 


1.4-3.1 


10.2 


8.2 


J/^ 




1.5 


1.3 


Hadrons 


3.1-12.0 


4.2 


2.8 


Hadrons 


> 12.0 


0.5 


0.2 


Total 




100.0 


100.0 



Table 5: Contributions to muonium HFS and a^ 



In Table ^ we compare results of various calculations of the leading order hadronic 
contribution to the muonium HFS. The first estimate of this effect was performed in p6| . 
The authors took into account the contributions of the p, uj and mesons and param- 
eterized the hadron continuum above 1 GeV under the assumption that R is constant 
and equals 2. This simplified approach gave nevertheless a result fairly close to those of 
the later, more sophisticated analyses. It is not surprising since the dominant part of the 
hadronic contribution comes just from the lowest vector mesons which properties were 



known quite accurately at that time. In |13 additional experimental data were taken 



into account above 1 GeV, including effects of heavy quarkonia. The hadron continuum 
was parameterized with a function R{s) = As^ in five energy ranges below 47 GeV 
and above this energy the asymptotic QCD formula with six quarks was used. Finally, 



in the recent work [|I^ the authors used a similar approach with a slightly more sophis- 
ticated parameterization of the hadron continuum in seven energy ranges below 60 GeV. 
Unfortunately, the authors of the above mentioned papers ignore the model dependence 
of their results which can be fairly strong and the error of their calculations only reflects 
the quality of the fit in their rather artificial models describing the data. Even larger is 
the effect of the systematic uncertainties completely ignored in most of the calculations. 

12 



Source Reference Ai^(had - vp), kHz 



J. R. Sapirstein et al, 1984 [0 0.22 ± 0.03 



A. Karimkhodzaev and R. N. Faustov, 1991 0.250 ± 0.016 



R. N. Faustov et al, 1999 [0 0.240 ± 0.007 

A. Czarnecki et al, 2001 This work 0.233 ± 0.003 

Table 6: Comparison of various calculations 



This is particularly true for the dominant contribution to the muonium HFS coming 
from the 27r channel. For example, in [|l^ its accuracy calculated within the model of [^ 



is unrealistically high and reaches 0.65% although the model itself is based on the data 
of ||66| where the systematic uncertainty in the dominant p meson region varies from 2 to 
4.4%. Direct integration of the experimental points in our work also assumes some model 
for the energy dependence of the data. We estimated the model dependence by comparing 
the trapezoidal integration in which experimental points are connected with straight lines 
with even simpler rectangular integration which assumes a constant cross section within 
a small energy range and found that both methods led to the same result, the difference 
being much smaller than the error. Both models allow a relatively simple estimation of 
the uncertainty arising because of the systematic errors. As already noted above, the 
effect is dominated by the contributions of the low lying vector mesons. Therefore, in all 
described calculations the central values of Ai/(had - vp) are close to each other. However, 
in our opinion their uncertainty might have been underestimated. The decrease of the 
uncertainty in the result presented in this paper became possible thanks to the utilization 
of the most recent data set coming from the high precision e~^e~ experiments in the low 
energy region. 

4 Conclusions 

The calculation of the hyperfine splitting in the ground state of muonium is summarized in 
Table |^. The biggest uncertainty of the theoretical expression comes from the inaccuracy 
of experimental determination of the muon mass and its magnetic moment. This inaccu- 
racy is essentially related to the statistical uncertainty of the present experiments 0, §, §] . 
Improvement of the statistical accuracy in muonium spectroscopy studies will be possible 
with future intense muon sources. Construction of such facilities is being considered in 
connection with muon storage rings which would serve as neutrino factories [^, ^. There 



are plans to build a high- intensity muon source at the Japanese Hadron Facility [|69[ . The 
muon production rates at those future sources are expected to be around 10^^ — 10^^ fi/s, 
higher than the most intense present beams (at the Paul Scherrer Institute) by a factor 
of 300-3000 or more. If muonium spectroscopy becomes part of the research program at 
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Correction Contribution to the HFS 

vf 4459031.920(511)(34) kHz 

Az/(ae) 5170.926(1) kHz 

Az/(QED3) -899.557 kHz 

Az/(QED4) -0.55(22) kHz 

Az/(weak) -0.065 kHz 

Azy(had - vp) 0.233(3) kHz 



Az/(had - h.o.) 


0.007(2) kHz 


^'HFs(theor) 


4 463 302.913(511) (34) (220) kHz 


i/HFs(exp) 


4463302.776(51) kHz 



Table 7: The muonium hyperfine structure: various theoretical contributions discussed 
in this paper 



these facilities, the statistical errors, which limit the present accuracy, will be decreased 



by one-two orders of magnitude |jTO[- In the more remote future, if a muon collider is 

constructed, further increase of beam intensity by a factor of about 100 will be possible. 

In the present paper we have studied the limits of possible bound-state QED tests with 

muonium hyperfine structure. Given the recent progress in multi-loop QED calculations 

2T, |7D|, one can expect that higher-order radiative effects will be evaluated when war- 



ranted by the experimental precision. However, the ultimate accuracy of every theoretical 
prediction will be limited by the knowledge of the hadronic effects. Here, we have studied 
the leading order vacuum polarization effects and found that they can be calculated with 
1.3% (or 3.1 Hz) accuracy. Higher order corrections can modify this contribution by a 
few percent. In the leading logarithmic approximation, they were found to increase the 



lowest order hadronic contribution by about 3% |T9|. In the case of the muon anoma- 
lous magnetic moment, this logarithmic correction is not dominant |7l|]. At present, it 
is not clear how good this analogy is and how reliable the logarithmic approximation 
is in the case of the muonium hyperfine splitting. When the experimental accuracy is 
improved, the higher order hadronic effects will have to be scrutinized. The evaluation 
of the diagrams with the vacuum polarization is straightforward, since only moderate 
relative accuracy is sufficient. The theoretical accuracy will be limited by the knowledge 
of the hadronic vacuum polarization in the leading order diagram, and by the ability to 
compute the hadronic light-by- light scattering diagrams which at present is possible only 
within models of low-energy hadronic interactions. 
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